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Abstract 

Semileptonic B decays into excited charmed mesons D' and D[ are studied in the framework of 
heavy quark effective field theory (HQEFT) up to order I/tuq. They are characterized by a single 
leading Isgur-Wise function r and several wave functions arising at 1/niQ order, r and the 1/niQ 
order functions Xq, Xq related to the kinetic energy operators are evaluated through QCD sum rule 
approach; zero recoil values of the 1/rriQ order functions n%, K2, k'i and k! 2 are extracted from the 
meson masses; and the branching ratios are found to be suppressed by the 1/mg corrections. It 
is concluded that the next leading order wave functions can be significant. However it does not 
change the previous prediction that the production rate of jf = | + charmed mesons dominates 
over that of ^ + doublets. 
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I. INTRODUCTION 



Semileptonic B decays are important in extracting the CKM matrix elements and ex- 
ploring CP violation. Presently the most promising approach to determine \V c b\ is to study 
either the inclusive semileptonic B decays or the exclusive decays to the ground state charmed 
mesons, B — > D{D*)lv. However the precision of these study depends on both the exper- 
imental measurements and the theoretical methods probing the nonperturbative effects of 
strong interaction. Semileptonic B decays into excited charmed mesons contain the main 
background for measuring the decays into D and D*, and they are important in relating the 
inclusive B decays to the sum of exclusive channels. To get precision knowledge on B physics 
it needs to study the decays into excited mesons from both experimental and theoretical 
aspects. 

For a hadron containing a single heavy quark Q (b or c) and any number of light quarks 
(u, d, s), the heavy quark spin sq is decoupled from the total angular momentum of the light 
degrees of freedom ji in the heavy quark limit ttiq — > oo. So ji becomes a good quantum 
number in this limit. Consequently, charmed mesons are usually classified by ji and parity. 
The ground state pseudoscalar and vector mesons (D, D*) have jf = \ . D x and D* 2 belong 

to the jf = | + doublet, while D' (or written as Dq in some references) and D[ are the | + 
one. 

In the recent years experiments made rapid progress on charmed meson spectroscopy, 
especially for the four lightest excited charmed mesons (Di, D 2 , D' and D[) as well as 
their counterparts of cs states. The | doublet mesons have narrow widths and their 
masses are known precisely: vnr> x — 2420MeV and mp* = 2460MeV Broad charmed 
resonances are observed in Dn and D*n systems by BELLE (2|, FOCUS [1J and CLEO 
Q Collaborations. The masses and widths of the | + doublets can be obtained from these 
measurements, nevertheless they still suffer from large uncertainties^. The branchingratios 
for semileptonic decays B -> D x {D* 2 )li> are reported by CLEO Q, ALEPH Q, DO @| and 
BELLE [9] Collaborations. Though not being confirmed, B — > D'^tv decay ratio is obtained 
recently by BELLE [9J. On the other hand these decays are studied by theorists via different 
approaches, among which are the operator product expansion (OPE) 0,113, Lattice QCD 
[l2l | and quark models [l3-15]. Note the I/tuq order corrections in the usual heavy quark 



expansion has been considered in the early work [14j, [15] . QCD sum rule method is also 
applied to calculate the form factors. Refs.[l6|, 17] studied the semileptonic B decays into 
excited charmed mesons at the leading order of heavy quark expansion (HQE), and the 1/tuq 
order contributions for B — > Di(D%)£v are considered in Refs. [l8| and [l9| using different 
framework of HQE. 

Generally, the theoretical calculations in the tuq — > 00 limit predict that the production 
of I doublets dominates over that of ^ doublets in semileptonic B decays. As illustrated 
in Ref. 11], for reasonable values of the Isgur-Wise function, the rate T(B — > D' (D[)£u) falls 
far below T(B — > D^D^iv). However, BELLE indicates in Ref.joj a large branching ratio 
for B decay to the wide D' Q state. If this result is confirmed, the previous throries or models 
need to be improved or corrected to explain it. Generally speaking, the predictions derived 
in the uiq — > 00 limit should always be supported by the estimation of 1/toq corrections, 
which turn out to be sizable in some specific situations. For example, it is known that the 



leptonic decay constants of heavy mesons receive considerable I/tuq corrections [20j, |21 
Calculations in different approaches also indicate large 1/mq corrections to the B — > D\k 



decay rate [15|, [18|, [19]. In the case of B — > D' (D' 1 )£i' transitions, one may ask whether the 
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great enhancement of the production rate for ^ + states is due to the finite mass corrections 
in the HQE. 

In this paper, the semileptonic B decays into the | + charmed meson doublet (D' , D[) 
are studied in the framework of HQEFT|22|-|25| that performs a complete decomposition of 
quantum fields and therefore includes the heavy quark- ant iquark coupling effects in the finite 
mass corrections. In SecJTT]we present the formulation of HQE to the decays B — > D'qID'^Iu. 
Up to the order of 1/rriQ, the relevant form factors are given by universal wave functions 
that are heavy flavor independent. In Sec lIIII QCD sum rule approach is used to evaluate 
the leading Isgur-Wise function £ and the next leading order wave functions Xo, Xo relevant 
to the kinetic energy operator. The QCD sum rule for leptonic decay constant of | + mesons 
is also derived. Sec JIVI analyzes the sum rules and gives numerical results. Finally a brief 
summary is given in SecJVl 



II. B D'^D'^Ip wave functions in hqeft 

The weak matrix elements relevant to B — > D' (D[)£u decays can be characterized by 
form factors as 

(D' (v')\cYb\B(v)} = 0, 
(D' (v')\c^ 5 b\B(v)) = ^/fn^mSig+i^ + v'n+g-i^-v'n), 
(DKv'^^crfblBiv)) = ^fn^mgv^ + igv^ + gvsv'n^-v), 
(D£(i/,e*)|c7V&|fl(v)) = ^ x m B %g A e^elv^' r (1) 

The initial and final states B and -Dq(I) are treated as heavy hadrons with the momentum 
■itlbv and mo' 0(1) v', respectively. The form factors gi are dimensionless functions of the 
product of velocities, y — v ■ v', and e* is the polarization vector of the axial vector meson 
D[. The differential decay rates are given by 

dT(B -+ D>£u) _ Gl\V cb \*m\ r r i{y2 _ l)3/2[(l + ^ _ (l _ ^ 



dy 487r 3 

" W - 1) 1/2 {2(1 - 2r[y + r?)[g 2 Vl + (y 2 - l)g\] 



dT(B^D[£u) G||K*| 2 m| m/3 , 2 _ , . r 2 , ,„2_ n „2 



dy 487T 3 

+[(y - r\)g Vl + (y 2 - l)(g Va + r[g V2 )} 2 } (3) 

m D i m D' 

with r' n = — 1 and r\ = — L . 

In the framework of HQEFT the matrix elements in QCD can be expanded in powers 
of 1/rriQ. Generally, the HQE of the matrix elements responsible for heavy meson leptonic 
decays and for transitions between heavy mesons can be written as [2c 



avy 
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^(0\qTQ\M) (0\qTQt\M v ) - ^<0|?T^(^ ± ) 2 Q+|M„> + 0(l/mj), (4) 



-(M'\Q'VQ\M) -> (M^\Q+TQ+\M V ) - — (M'JQ+T—(iI? ± ) 2 Q+\M v ) 



- J_(M;|g+(-^) 2 ^rg+|M v ) + o(i/m 2 Q ) (5) 



with the definition 



Djj 4 = v ^v ■ D, (6) 
D^ = D^~ v»v ■ D, (7) 

K^ip = - I KD^ip. (8) 



M(M') can be any ground or excited heavy meson containing a single heavy quark. Q is 
the field in the full QCD Lagrangian, and is the effective heavy quark field in HQEFT, 
carrying only the residual momentum k = pq — mqv. \M) is the meson state in the full 
theory, while \M V ) is an effective state defined in the HQEFT so as to display the heavy 
quark spin-flavor symmetry. They are normalized as 

(M\QYQ\M) = 2m M v^ (9) 
(M v \Q+>fQ+\M v ) = 2Av", (10) 

where Am = thm — wq is the mass difference between the heavy meson and heavy quark, 
and 

A = lim Ah = lim (mjy ~ m o) (11) 

mQ->oo toq— >oo 

is the heavy flavor independent binding energy. The state \M V ) defined in this way is 
irrelevant to the heavy quark mass and related to \M) via 

(M'\Q'TQ\M) = 1 _ {M' vl \J eff e^ dixC ^i\M v ), (12) 



y/m M ,m M \fA M ,A 



M 



where £ e // is the HQEFT Lagrangian and J e // is the effective current for QTQ, gene rally 



also written as an expansion in I/tuq. For more details of HQEFT we refer to Refs. |22H25 

Due to the heavy quark symmetry, form factors for heavy-to-heavy transition matrix 
elements can be parameterized by a set of wave functions, which are universal in that they 



are heavy flavor and spin independent. To define these wave functions one may use the 

'i ) and excited (| 



following spin wave functions for the ground (\ ) and excited (| ) states: 



k. = va>p + \ _^ { Z $| (M) 

where A = Ai- and A' = Ai+ are the binding energies of the h and ^ + doublets, respec- 

2 2 ' Z Z 

tively, and P± = (1± yf^)/2. Then in HQEFT the matrix elements between | + states can 
be parameterized as 

(K v ,\Q+j>*Qt\K v ) = e(y)Tr[jC v/ j»lC v }, 
(K v ,\Q^j^~a a pF a PQt\K v ) = ^Tr^^v'^rP^-a^] (15) 
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with K v i = 7°/C^,7°. The states K v and K v > can be either of the two mesons belonging to 
the | + doublet. The tensor n' al3 (v,v') are decomposed as 

Kf°P{v, v') = m' 2 a a/3 - 4(7 Q v ,/3 - 7V ) + k'^v? - 7 V) + K,' 5 (v a v'P - v p v' a ) (16) 

with £' and k[ being scalar functions of y = v ■ v'. Carrying out the trace calculation in ( fl5|) 
and setting v' = v, one gets from Eqs.©, <Q, ( ITOj) . (fT2l) and (IT51) 

jjm^y = ^{2A'e'(i) - — K(l) + 34(1))K, (17) 
W y = ^i{_ 2 A'£'(l) + — K(l) - 4(l))}( e * ■ eK. (18) 

Since the Isgur-Wise function £'(1) satisfies the normalization condition £'(1) = 1; * ne above 
equations yield 

Ajy o = A'- — («i(l) + 3«4(1)), (19) 

77l c 

A Di = A' —(k[(1) — «4(1))> (20) 

77l c 

which are quite similar to those relations for the jf = | ground state bottom mesons 



23|,|24 
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Ab = A-— (ki(1) + 3k 2 (1)), (21) 
m b 

Ab* = A-— ( Ki (1)-k 2 (1)). (22) 

For I? — > D' (D[) transitions, the relevant matrix elements can be parameterized in 
HQEFT as 

{K V ,\Q+TQ+\B V ) = T{y)Tr[£ v ,TM v l 
(K V ,\Q+T-^D 2 ± Q+\B V ) = - X h (y)jTr[iC v ,TM v ], 

(K v ,\Q+Dl—^TQt\B v ) = - x c( y )l-Tr[iC v ,rM v }, 
-iv'- D A 

(K v ,\Q+T-^ l -a a pF a PQ+\B v ) = -lTr[R b a p(v,v')iC v ,rP + ia a pM v ], 

{K v ,\Q+ l -a a pF a P-^TQ+\B v ) = -l-Tr[R^(v,v')jC v >ia^P[.rM v }, (23) 
1 -iv'- D A 

where the Lorentz tensors R b ^\v,v') can be decomposed as 

Btp{v, v') = xtlajp + xKjp, (24) 
Kp( v i v ') = xllcTIa + X C 2V a lfi- (25) 
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The wave functions r and Xi (i = 0, 1, 2) depend on y. r is dimensionless, while Xi nas 
mass dimension two. 

HQE for the form factors Qi is then obtained from ([1]), ([5]) and f l23|) . Up to the order of 
1/mQ one has 

g+ = o, 

1 7 [ 3 Xi-xS(i + y)], 



with 



m c A 

9Vl = [f + 2^ (Kl(1) + 3K2(1)) + 2^ K(1)-ft:,2(1)) "i xb 

1 ; x!](y-i), 



fiV 5 



m c A 

X% 



m r A 



m c A' 

9 - = f + 2i (K ' (1) + 3K2(1)) + 2^ W(1) " ~ ^ 

1 ~,X\ (26) 



m r .A 



b c 

f = T (27) 

2m b A 2m c A'' v y 

x b = 3x b i-(i + y)xl (28) 

Here /cf } (l) (i = 1,2) are the zero recoil values of , whereas other wave functions and 
form factors depend on the variable y = v ■ v '. 



III. QCD SUM RULES FOR WAVE FUNCTIONS 

As can be seen in ( )26|) . in the heavy quark limit all form factors simply reduce to the 
Isgur-Wise function r. Among the 6 functions x^ C \^ = 0,1,2) of order I/tuq, Xi^) are 
defined in ( )23|) by the chromomagnetic operators. Contributions from such operators are 



generally expected to be v ery small, which is supported by the relativistic quark model [26 
and QCD sum rule study [27|. Here we mainly focus on the functions Xo an d Xo> which are 
defined by the matrix elements of the kinetic energy operators. Since the kinetic operators 
preserve heavy quark spin symmetry, Xo simply correct the leading Isgur-Wise function r 
in the way of Eq. ([2"7]) . 
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In order to calculate r and Xo> Xo> we study the following three-point correlation functions 
S T = i 2 I dWze^'^OlTl^^ (29) 



~*o = e J *arfSe^*^ (30) 

/ d 4 xd 4 ^( fc '-^(0|T{J 0)+ , 1/2 (x), (Q+ Sf— ^-rO+)(0), jj_ 1/2 (*)}|0>(31) 

-zV- D 

where k and fc' are the residual momenta of the heavy quarks. T should be 7 M and 7^7 5 for 
vector and axial vector heavy quark currents respectively. Jj t pj t with j the total spin of the 
meson should be proper interpolating currents for the heavy-light mesons. One set of such 
currents are proposed in Ref.|28j. One can use 

4,-,i/2 = ^Qt^Q, (32) 

4,-,i/2 = \fl0hlQ, (33) 
for the | ground state doublet, and 

4, + ,i/2 = \I\qU, (34) 

4, + ,i/2 = (35) 

or 

Jl + M2 = \J\Q + A-^^ (36) 
rt _ . / I/o+„,5„,a, 



Jl, +l i/2=\I^QZ'rilHW±Q (37) 

for the | + doublet. 7" is defined as 7" = 7° — v a ft. In Eqs. (l29l) - (I3~T]) Jq,+,i/2 is used in 
the three-point functions. Of course one can substitute Ji,+,i/2 for Jo,+,i/2 in the evaluation, 

which does not make difference to the results for wave functions r and Xo ■> as required by 
the heavy quark symmetry. 

The formulae in (129 p - ( 13 ip are analytic functions of the variables u = 2v-k and u>' = 2v'-k' 
with discontinuities for their positive values. The phenomenological representation for these 
correlators can be obtained by inserting the complete set of intermediate states with the 
same quantum numbers as the currents Jo,+,i/2 and Jo,-, 1/2- Isolating the pole terms of the 
lowest states we get 

(0|J 0f+4/a |l^>(l^|Q}rQ+|S t ,>(B 1 ,|jJ_ 1/!J |0> 
^phen = 7 — T x" ; \ a a) ^ n x S ner resonances, (38) 



(2A - w - z'e)(2A' - w' - z'e)AA' 
(2A - cu - z'e)(2A' - u' - z'e)AA' 



"phen = 777a 77 — 7777777T7 777 — „• A x a / ^ ill gJ ier resonances, (39) 



(0 1 Jo,+,i/2 1 K> ) I Dl-^TQ+ \B V )(B V \ J Q T _ 1/2 1 0) 
= (2A-c-z e )(2A"-^-z e )AA + higher resonance440) 



where the first term in each equation is a double-pole contribution, and the second term 
takes into account higher states and continuum contributions. Using T = -y^ 5 and the 
definition in ([23"]) . one gets the pole terms: 



' — 'pole 



(2A - u - ie)(2A' - u' - ie) 



v- v y, (41) 

v - v'f, (42) 



-pole ~ (2A-u- ie )(2A' - a/ - ie) A 

£(«-«T- (43) 



^pole 



(2A - w - ze)(2A' - fa/ - ie) A' 



/i+ and /i- are the leptonic decay constants of relevant mesons at leading order approxi- 

2 2 

mat ion: 



(0|J , + ,i/ 2 |£>o,/> = VA'J+,1/2, (44) 
(0\J ^, 1/2 \B v ) = v/A/Li/2- (45) 

In sum rule approach, the theoretical representation for the correlation functions can be 
calculated from QCD or effective theories in the deep Euclidean region, and in performing the 
operator product expansion the nonperturbative effects are incorporated via the inclusion 
of nonzero vacuum condensate values. Formally the theoretical sides of the sum rules can 
be written as 



6(c) r p-r t (p b ^ t ) 

-//„ o\^ X tLo ) = J dvdv ' -r - - r^TZ u> - it) + Enp + subtraction terms ( 46 ) 



with E NP being the nonperturbative terms. QCD sum rules are obtained by equating 
the phenomenological and theoretical representations of the correlators. In doing this the 
perturbative contribution above some threshold energy is assumed to simulate the higher 
resonance contribution. To suppress the higher resonance contribution and at the same time 
enhance the importance of low dimension condensates, Borel transformation 

BP = T lim -£_(_- 1)" with T = — fixed (47) 
n-too-oj-^-oo i{n) du n 

should be performed to both sides of sum rules. Since there are two variables u and u' for the 
correlation functions (|29l) -( l3Tl) . we shall perform a double Borel transformation B^B^,, 
which then introduces two Borel parameters t and t' in the sum rules. In studying B decays 
into ground state charmed mesons, it is argued that the hadronic and perturbative 



spectral densities can not be locally dual to each other, but the quark-hadron duality is 
restored in the "diagonal" variable v + = ^y l - Here we shall follow this prescription. That 
is, we integrate the spectral densities over the "off-diagonal" variable z/_ = and assume 
the quark-hadron duality in u + for the integrated spectral densities. This can be represented 
as 

Zpoie = 2 / dv + e~ u+/T p pert (is + ) + E NP , (48) 
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where the two Borel parameters are set equal, t = t' = 2T. H is obtained by applying double 
Borel operators to H, and 



PpertKy^t 



dv_p pert (v + ,v-). (49) 



In the OPE we consider condensates with dimension no higher than 5, and the light quark 
mass and higher radiative corrections are neglected. Then the Feynman diagrams presented 
in Fig.l should be calculated. The resulting sum rules turn out to be 

/i + A-re- (A+A ' )/T = , 2 * V2 C du + ule-^- 2 ^{q q ) + -^\^{y + 2) 

2 2 8ll 2 {l + y) 2 J 671 9bTT 2 l 



37r(y + l) 



f f Xo -(A+A')/T _ y + 4 f S ° , A -u+/T 

a 6 - 16^(1+^ y dv+v + e 

+ 9 67 r(l + ^^ (F ^ F >' 
(y + 8)T 



(51) 



^a s (F*F aa P} } (52) 
967r(y + l) sV «/3 /' v ; 

where the threshold values should be determined by the principle of minimal sensitivity in 
the numerical analysis of sum rules. The condensates have the typical values: 

(qq) w -(0.23 GeV) 3 , 

%(qa afi F afi q) « -m 2 (qq) with m 2 = 0.8 GeV 2 , 

a s (F^F aaf} ) « 0.04 GeV 4 . (53) 

Eqs. ( 13"2~j) and ( 13T)j) are used as interpolating currents in deriving the sum rules (150]) - (152"]) . 
We have also considered the current (I34p but we find that using such current results in 
zero contribution of the perturbative diagram (the first diagram in Fig.l), which makes the 
resulting sum rule equation not reliable. This has been noted in Ref.[l7[ and we just further 
checked it. 

To derive the wave functions from these sum rules one needs to know the leptonic decay 
constants f\- and fi + . They can also be evaluated in the same framework through QCD 

sum rule approach. The sum rule for f L - has been analyzed by previous work 20, [24], 3l| 

and our result is [2 



2 l07T z J I 67T 8i z 7T 



48vrT 



a 



S (F^F^), (54) 



where the relation between f L - and F in Ref . 24] is F = y/2f L - . 
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For + , we consider the two-point correlation function 

2 



n = i J d 4 xe lk - x (0\T{J ^ 1/2 (x), Jj +jl/2 (0)}|0). (55) 

Inserting a complete set of intermediate states and assuming the quark-hadron duality, one 
has 



fl 



dv ^ pert ; — h H-NP + subtraction terms. (56) 



2A' — 2v ■ k — ie J v — u> — it 
After Borel transformation we get 

/2 -2A'/T = 3 f'o 4 U/T ( _3_ _ Os_^ F aP y (57) 

s 64tt 2 7 v 16 32vr w v ; 

/i+ has also been studied in the usual HQET j28[. We note that the perturbation term in 

Eq. floTl) is same as that in Ref.[28j]. Our calculation includes contributions from all diagrams 
in Fig. 2, and the nonperturbative terms have some difference to that reference. 

Sum rules in Eqs. fl5T?]) - fl52l . f[54"j) and (1571) constitute the main results that we will use to 

discuss the B — > D'^D'^iv decays. The constants f k - and fi+ as well as the binding energy 

_ _ 22 

A and A' can be estimated from sum rules (1511) and f|5T|) . And the wave functions r(y), 
Xo(y) and Xoiv) can derived by studying the ratios of Eqs. (l5T!|) - (l52"|) to (1511) and fl5T|) . 

QCD higher order corrections are not included in our calculation. They affect both the 
three-point and two-point correlation functions and deserve further study in future work. As 
far as the determination of transition wave functions is concerned in this paper, the effects 
of radiative corrections are expected to be partly cancelled in the ratios of three-point to 
two-point correlators, and therefore not influence the final results significantly. This has 



been proved to be true in the study of Refs. [16|, |32j, |33| . In those references the two- loop 



corrections to the Isgur-Wise functions are found to be small and well under control for the 



B decays into both ground state [32|, |33| and excited state [16j charmed mesons, although 



the corrections to decay constants are sizable. 

IV. NUMERICAL RESULTS 

We get from Eq. (1541) the appropriate binding energy and decay constant as 24 

Ai- = 0.53 ± 0.08GeV, 

2 

/ x - = 0.21 ±0.05GeV 3/2 . (58) 

2 

For fi+ one should study Eq. (l57|) . Ai+ and f\+ as functions of the Borel parameter 

2 _ 2 2 

T is presented in Fig. 3. Ai+ and f\+ have acceptable stability when setting the threshold 

2 2 

Sq = 2.6 — 3.0GeV. The curves in Fig. 3 become rather stable when T > lGeV. However, the 
criterion of sum rule analysis is that both contributions from the higher resonances and from 
the higher order power corrections in OPE should not be very large, say not much higher than 
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30%. According to this criterion the proper window for Eq. (l5"T|) is 0.6GeV < T < 0.8GeV. 
As a result we get 

Ai+ = 0.81 ±0.12GeV, 

2 

f 1+ = 0.30 ± 0.05GeV 5/2 , (59) 

2 

where the central values are obtained using s$ = 2.8GeV and T = 0.7GeV, and the errors 
are attributed to the variation of the threshold and Borel parameter. 

Now we can study the sum rules for B — > D' (D[)£u wave functions. The leading function 
r depends on the recoil variable y and can be estimated from the sum rule (l50l . Fig. 4 displays 
r as a function of the Borel parameter at the fixed point y = 1. Applying the sum rule 
criterion the appropriate region for analyzing the stability is 0.8GeV < T < l.OGeV. As 
can be seen in the figure, 2.7GeV < Sq < 3.3GeV is favorable. Therefore we get 

r(l) = 0.57 ±0.06, (60) 

where the central value corresponds to Sq = 3.0GeV and T = 0.9GeV. 

Following the same procedure the subleading order wave functions Xo an d Xo can be 
derived from Eqs. fl5T|) and fl52l) . The results at zero recoil are shown in Fig. 5 and Fig. 6, 
respectively. In the appropriate window Xo(l) * s n °t sensitive to the Borel parameter when 
Sq ~ 2.1GeV, while Xo(l) becomes stable around a smaller threshold value Sq ~ 1.2GeV. 
Setting T ~ 0.9GeV we then obtain the following zero recoil values for the 1/rriQ order wave 
functions 



xS(i) 

A 

xS(i) 

A' 



0.27±0.12GeV, (61) 
-0.20±0.12GeV, (62) 



where the errors mainly arise from the thresholds. So Xo an d Xo have opposite signs. As 
nib > m c , Xo ma y yield a relatively larger contribution to the B — > D' Q (D[) form factors. 
Xo can only weakly counteract the contribution of Xo, which makes f in Eq.( j27|) suppressed 
with respect to r. 

If one fix the values of the thresholds and the parameter T, r and Xo as functions of the 
recoil value can be evaluated from the sum rule equations. The results are shown in Fig. 7, 
where T = 0.9GeV is used. We find these functions can be expanded near y — 1 as 

r(y) = r(l)[l - 0.56(y - 1) + 0.35(y - l) 2 ], (63) 
X b (y) = Xo(l)[l - 1-45(1/ - 1) + 0.98(1/ - l) 2 ], (64) 
Xo(v) = Xo(l)[l " 0.52(y - 1) + 0.28(y - l) 2 ]. (65) 

When Xi(2) are neglected, the form factors in (T2"6"|) can be simply written as 

9+ = 0, g- = Tjy o , g Vl = (y- l)r D > v 

9V 2 = 0, g Vi = -f D , v g A = T D i i (66) 

with the definition 

** =f + 5^A (Kl(1) + 3K2(1)) + 2^ W(1) + 3K ' M) - (67) 
= f + ^A^ 1 * 1 ' + 3K2(1)) + 2^> (<(l) - 4<1)) ' <68) 
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Consequently the differential decay rates turn into 



dT(B -» D' Q £i?) G 2 F \V cb \ 2 m% 3 2 , 3/2 2 

^ " 487T3 ^ (1 - ^ ~ lj TD 'o ' 

^_^) = ^W ri 3 (y2 _ 1)1/2[(1 + r;2)(5 , 2 _ %+1) 

-rU% 3 -10y 2 ±4</-2)]f 2 ,. (69) 

«f (l)(i = 1,2) are parameters related to meson masses. Taking m& = 4.67 ± 0.05GeV, 
m c = 1.35 ± 0.05GeV, rriB = 5.279GeV, mg. = 5.325GeV and the averaged masses for 
the \ + doublet, ra D / = 2.351 ± 0.027GeV and m D[ = 2.438 ± 0.030GeV 0, we get from 
Eqs.p}-fl22P 

Kl (l) = -0.53 ± 0.23GeV 2 , k 2 (1) = 0.05 ± O.OlGeV 2 , 

^(1) = -0.35±0.17GeV 2 , 4(1) = 0.03 ± O.OlGeV 2 , (70) 

where we include uncertainties from the binding energies and the quark and meson masses. 
The values of «i(l) and k 2 (1) in (I70|) are consistent with the results of Ref.(24]| in which 
these two parameters are evaluated through QCD sum rule equations that are independent 
of the heavy quark and meson masses. k\ (1) characterize the mass splittings of two mesons 
belonging to a jf doublet, and the absolute values of them are much smaller than those of 



K 



1). Therefore one gets from Eqs.()67l) and ( |68l) tjji ~ f D > L . It is also clear from ( !70l) that 

Try and r£>/ may be further suppressed by Krp with respect to f and r. 

With the obtained values of wave functions, we get the decay rates and branching ratios 
in Table 1. In the calculation the B meson life time t(B) = 1.542 ps and \V c b\ = 0.041 are 
used, and the uncertainties arise from the quark and meson masses as well as the variation 
of the thresholds and the Borel parameter. It is shown that the branching ratios in the 
itlq — > oo limit can be about 1 x 10~ 3 . However Xo an d Xo ma Y lead to nearly 20% 
suppression. When the contribution of «j (1) is included, the branching ratios can even 
be significantly reduced. We note that the masses of | + charmed doublets have not been 
determined precisely. Consequently the values of may suffer from larger uncertainty. 
So do the data in the last column of Table 1. Nevertheless the substantial suppression effect 
of the 1/rriQ contribution to the decay rates is evident. 







rriQ — > oo limit 


with 1/rriQ correction 
from Xo 


with 1/m.Q correction 
from Xo an d K T 


B D' £u 


r 


4.14 ± 1.20 


3.45 ± 1.09 


2.13 ±0.69 


Br 


0.97 ±0.28 


0.81 ±0.26 


0.50 ±0.16 


B D[iv 


r 


4.65 ± 1.32 


3.89 ± 1.20 


2.05 ±0.66 


Br 


1.09 ±0.31 


0.92 ±0.28 


0.48 ±0.16 



Table 1. Rates T (in units of |V cb /0.041| 2 x 1CT 16 GeV) and branching 
ratios (in 10~ 3 ) of B — > D' (D' 1 )£p decays in the tuq — > oo limit as well 

as when taking account of 1/mg order corrections from Xo an d K i\^)- 
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V. SUMMARY 



Semileptonic B decays into jf = | + doublet excited charmed mesons are studied in the 
framework of heavy quark effective field theory with inclusion of the heavy quark- ant iquark 
coupling effects in the finite mass corrections. We present the heavy quark expansion for B — > 
D' (D[) transition matrix elements. At the leading order of HQE all form factors reduce to 

the leading Isgur-Wise function r, and at 1/ vtlq order there are six wave functions (i = 
0, 1,2). Among them Xo characterize the contribution from the kinetic energy operators 
and are expected to be much larger than Xi(2) that are relied to the chromomagnetic 
operators. Beside these functions characterizing transitions between | and \ + doublets, 

we present also the functions k[{i = 1, 5) for transitions between ^ + mesons. The zero recoil 
values of k[ and extracted from the excited meson masses. 

QCD sum rule method is applied to evaluate the functions r, Xo an d Xo> an d the decay 
rates are predicted. The Isgur-Wise function r gives the branching ratios of the magnitude 
10 -3 . However, the ratios are suppressed rather than enhanced by the 1/toq corrections. 
Though we have not calculated all 1/mg order wave functions, our results imply that the 
finite mass corrections would not likely change the dominance of the semileptonic B decay 
rates to | + states over the rates to \ + states. As a result, the production of | (D,D*), 

| + (L>i,L>2) and | + (-Dq,.D^) charmed mesons do not saturate the total semileptonic decay 
rate Fsl{B), and the configuration for the "missing rate" remains an interesting question. 
This is in agreement with the conclusion of the early work 15[ that adopted some model 
dependent assumptions. 

As for the recent report of BELLE [q| there is no indication of a broad D[ in the 
B — > D*tt£u channel, but the measurements indicate for B — > D' £u a large rate of similar 
magnitude to the § production rates. If that result is confirmed, the framework to predict 

| + production should be improved to connect the theories and measurements. Anyway, 
according to the calculations those rates might be in the reach of B facilities. Measurements 
on such processes will test the theories and shed light on the nature of excited states. 



Acknowledgments 

The author would like to thank Prof. Y. L. Wu for stimulating discussions. This work 
was supported in part by the National Science Foundation of China (NSFC) under the grant 
No. 10805005; and by the Project of Knowledge Innovation Program (PKIP) of Chinese 
Academy of Sciences, Grant No. KJCX2.YW.W10. 



[1] Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 1 (2008). 

[2] BELLE Collaboration, K. Abe et al, Phys. Rev. D 69, 112002 (2004). 

[3] FOCUS Collaboration, J. M. Link et al, Phys. Lett. B 586, 11 (2004). 

[4] CLEO Collaboration, S. Anderson et al, Nucl. Phys. A 663, 647 (2000). 

[5] P. Colangelo, F. De Fazio and R. Ferrandes, Mod. Phys. Lett. A 19 2083 (2004). 

[6] CLEO Collaboration, A. Anastassov et al, Phys. Rev. Lett. 80, 4127 (1998). 

[7] ALEPH Collaboration, D. Buskulic et al, Z. Phys. C 73, 601 (1997). 



13 



[8] DO Collaboration, V. M. Abazov et al, Phys. Rev. Lett. 95, 171803 (2005). 

[9] BELLE Collaboration, D. Liventsev et al, Phys. Rev. D 77, 091503 (2008). 

[10] N. Uraltsev, Phys. Lett. B 501, 86 (2001). 

[11] I. I. Bigi, B. Blossier, A. L. Yaouanc, L. Oliver, O. Pene, J. C. Raynal, A. Oyanguren and P. 

Roudeau, Eur. Phys. J. C 52, 975 (2007). 

[12] D. Becirevic et al, Phys. Lett. B 609, 298 (2005); B. Blossier, A. L. Yaouanc, V. Morenas 

and O. Pene, Phys. Lett. B 632, 319 (2006); Erratum ibid 645, 476 (2007). 

[13] F. Jugeau, A. L. Yaouanc, L. Oliver and J. C. Raynal, Phys. Rev. D 72, 094010 (2005). 

[14] A. K. Lcibovich, Z. Ligeti, I. W. Stewart and M. B. Wise, Phys. Rev. Lett. 78 3995 (1997). 

[15] A. K. Leibovich, Z. Ligeti, I. W. Stewart and M. B. Wise, Phys. Rev. D 57, 308 (1998). 

[16] P. Colangelo, F. D. Fazio and N. Paver, Phys. Rev. D 58, 116005 (1998). 

[17] Y. B. Dai and M. Q. Huang, Phys. Rev. D 59, 034018 (1999). 

[18] M. Q. Huang and Y. B. Dai, Phys. Rev. D 64, 014034 (2001). 

[19] W. Y. Wang and Y. L. Wu, Int. J. Mod. Phys. A 16, 2505 (2001). 

[20] M. Neubert, Phys. Rev. D 46, 1076 (1992). 

[21] J. M. Flynn and C. T. Sachrajda, Adv. Ser. Direct. High Energy Phys. 15, 402 (1998). 

[22] Y. L. Wu, Mod. Phys. Lett. A 8, 819 (1993). 

[23] W. Y. Wang, Y. L. Wu and Y. A. Yan, Int. J. Mod. Phys. A 15, 1817 (2000). 

[24] W. Y. Wang and Y. L. Wu, Int. J. Mod. Phys. A 16, 377 (2001). 

[25] Y. L. Wu, Int. J. Mod. Phys. A 21, 5743 (2006). 

[26] D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D 61, 014016 (2000); 62, 014032 
(2000). 

[27] M. Neubert, Z. Ligeti and Y. Nir, Phys. Lett. B 301, 101 (1993); Phys. Rev. D 47, 5060 
(1993). 

[28] Y. B. Dai, C. S. Huang, M. Q. Huang and C. Liu, Phys. Lett. B 390, 350 (1997). 

[29] M. Neubert, Phys. Rev. D 46, 3914 (1992). 

[30] B. Blok and M. Shifman, Phys. Rev. D 47, 2949 (1993). 

[31] M. Neubert, Phys. Rep. 245, 259 (1994). 

[32] M. Neubert, Phys. Rev. D 47, 4063 (1993). 

[33] E. Bagan, P. Ball and P. Gosdzinsky, Phys. Lett. B 301, 249 (1993). 



14 




Fig.l. Feynman diagrams contributing to the sum rules for r, Xq an d X(j- 
The thick lines represent heavy quarks; the light lines are light quarks; 
the curves are gluon fields; the black dots represent condensates; and 
the external lines represent the currents in Eqs. (l29p - (l3ip . 
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Fig. 2. Feynman diagrams contributing to the sum rule for + . The 

2 

external dashed lines represent the interpolating currents in Eq. (|55p . 
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Fig.3. Ai+ (figure (a)) and /i+ (figure (b)) as functions of Borel pa- 

2 2 

rameter T. The dashed, solid and dotted curves correspond to s^=2.6, 
2.8 and 3.0 GeV, respectively. 




Fig. 4. t(1) as a function of the Borel parameter T. The dashed, solid 
and dotted curves correspond to Sq=2.7, 3.0 and 3.3 GeV, respectively. 
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Fig.5. — Xq(1)/A as a function of the Borel parameter T. The dashed, 
solid and dotted curves correspond to Sq=1.9, 2.1 and 2.3 GeV, respec- 
tively. 



0.2 

D U 



- -0.2 




~-0.4 ! > , " : 

H 

u £-0.6- 
i 

-0.8- 

0.25 0.5 0.75 1 1.25 1.5 1.75 

T (GeV) 

Fig. 6. — Xq(1)/A' as a function of the Borel parameter T. The dashed, 
solid and dotted curves correspond to Sq=0.9, 1.2 and 1.4 GeV, respec- 
tively. 
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Fig. 7. r, — Xo/A and — Xo/A' as functions of the variable y. The dashed, 
solid and dotted curves correspond to Sq=2.7, 3.0 and 3.3 GeV in (a); 
s b =1.9, 2.1 and 2.3 GeV in (b); and sg=0.9, 1.2 and 1.4 GeV in (c). In 
the evaluation the Borel parameter T is set to 0.9 GeV. 
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